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PREFACE 

In  this  note  we  describe  and  collect  for  future  reference 
in  this  series  some  general  properties  of  the  equations  of  fluid 
magnetics.   The  generality  of  the  equations  described  in  MH-1  makes 
mathematical  analysis  forbidding,  and  it  is  necessary  to  specialize 
for  the  purpose  of  analysis.   In  this  supplement  to  MH-I  we  describe 
the  fluid  magnetic  systems  to  which  we  shall  have  most  frequent 
reference.   Here,  we  shall  adopt  a  set  of  standard  non-dissipative 
fluid  equations  coupled  with  Maxwell's  equations  by  a  Lorentz  force 
term  and  some  form  of  Ohm's  law.   Within  this  frame,  there  is  a  con- 
siderable amoiont  of  flexibility;  e.g.,  in  the  choice  of  an  equation 
of  state  and  of  a  form  of  Ohm's  law,  and  in  the  decision  to  include 
or  omit  displacement  current,  electrostatic  forces,  and  electro- 
magnetic momentum.   An  important  criterion  in  making  these  decisions 
is  the  mathematical  tractability  of  the  resultant  equations;  this 
desideratum  should  on  occasion  even  override  in  importance  the 

criterion  of  the  precise  physical  validity  of  a  given  assumption. 

The  assistance  of  H.  Kranzer  in  preparing  this  note  is  gratefully 

acknowledged , 
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1,   The  Electromagnetic  Equations 

To  describe  the  electromagnetic  field  we  shall  employ 
either  the  complete  equations  of  Maxwell 


'   15  +  curl  E  =  0,   div  B  =  0 
dt  ' 


(1.1)  } 

-  —  II  +  curl  B  =  [iJ,      div  E  =  qA 
c 

or,  alternatively,  the  system  in  which  displacement  current  is 
omitted, 

r  3B 

(1.2) 


dt 


+   curl  E  =  0,    div  B  =  0 


I  curl  B  =  p-J,   div  E  =  q/K 


This  system  vjill  be  called  the  pre-Maxwell  equations  since  it  is 
a  formalization  of  the  electromagnetic  theory  in  existence  before 
Maxwell's  introduction  of  displacement  current.   Maxwell's 
equations  are  Lorentz  invariant,  but  the  system  (1.2)  is 
invariant  under  Galilean  transformation  with  the  quantities 
E  =  E  +  uvB  (polar)  and  B  (axial)  treated  as  vectors.   Since 
we  shall  wish  to  couple  the  electromagnetic  equations  v;ith  non- 
relativistic  equations  of  the  fluid, the  system  (1.2)  is,  in  a 
sense,  the  more  appropriate,  but  for  the  present  the  choice  of 
electromagnetic  eouations  is  left  open.   Prom  (1.1)  we  obtain 
the  conventional  form  of  the  charge  conservation  equation 

See  MH-V,  Theory  of  Maxwell's  Equations  without  Displacement 
Current. 
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(1.3)    1^  +  div  J  =  0 

but  from  (1,2)  we  have 

(l.li)     dlv  J  =  0. 

Accordingly,  the  rystein  (1.2)  admits   only  solutions  correspond- 
ing to  the  flovj  of  currents  in  closed  circuits;  in  particular, 
the  singular  case  of  a  moving  point  charge  cannot  be  considered 
as  a  localized  singular  current  source. 

The  energy  equation  derived  from  (1.1)  is 

Prom  (1.2)  we  obtain  instead  of  (1.5)  the  relation 


^  B^)  +  div(i 


(1.6)     ^  (-^  B^)  +  div(i  E/B)  =  -E.J, 


and  this  will  be  taJcen  as  the  formulation  of  energy  conservation 
for  the  system.   In  other  words,  v.'hen  we  drop  displacement 
current  we  treat  magnetic  energy  alone;  the  significance  of  this 
will  later  be  made  clear. 

The  electromagnetic  momentum  equation  may  be  obtained  from 
(1.1)  by  Use  of  the  identity 

2 
(1.7)  curl   ZX  Z  =  div(Z^Z^   -  -^r  S^"^')    +   Z^div  Z 


as  follows: 


(1,6)  -x4r(EXB)    +    dlv   M  =   qE  +    J^B 
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where  the  Maxwell  stress  tensor  is  defined  by 


M  =  M^  +  I^ 


uij  =x\eV    -   i  E^S^J 


(1.9)  <  M^'  =/^\E"E^  -  I  E^5- 

Using  the  identity  (1.7)  vje  obtain  from  the  system  (1.2)  the 
equation 

(1.10)  div  Mg   =   J  X  B 

which  ^^;ill  be  taken  as  the  appropriate  formulation  of  momentum 
conservation  for  the  pre-Maxwell  system.   V/hen  we  ignore  displace- 
ment current,  then,  we  choose  to  ignore  electric  forces  as  compared 
with  magnetic  forces  as  well  as  electric  energy  compared  with 
magnetic  energy.   These  stipulations  will  appear  less  arbitrary  and 
even  necessary  when  we  consider  the  total  energy  and  momentum  equa- 
tions including  both  the  fluid  and  the  electromagnetic  contributions, 


Integral  form  of  the  electromagnetic  equations 

We  shall  treat  problems  for  which  there  may  occur  singular  sur- 
faces where  we  allow  discontinuities  in  certain  quantities.   We  shall 
therefore  take  the  int'^gral  form  of  the  equations  as  f\andamental  and 
derive  discontinuity  conditions  from  these.-   Maxwell's  equations  may 
be  written  in  integral  form  as  follows: 
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(1.11a)   ^    B.dS 
Z 


E'dx  =  0, 


(1.11b) 


B-dS  =   0, 


(1.11c)    -  ^ 


1_  _d_ 
2  dt 


E«dS  +  q)  B'dx  =  [i 

r 


J'dS, 


J, 


(l.lld)   C)  E'dS 

s 


1  f 


q  dV, 


D 


In  these  equations  P  stands  for  the  boundary  curve  of  a  surface 
y~,    S  for  the  boundary  surface  of  a  domain  D.   Throughout,  \^e 
shall  use  the  symbols  D  to  denote  a  domain,  S  to  denote  a  closed 
surface,  ^I  ^  surface  with  boundary,  P  a  closed  curve,  G  an  open 
arc.   For  moving  surfaces  we  have  the  alternative  forms. 


r 


(1-12)   ^ 


B-dS   + 

'zct) 


)    (E  +  v  X  B)  •  dx   =   0 
P(t) 


(1.13) 


1   d 
c 


E-dS  + 

£(t) 


>    (B 

r(t) 


c 


V  /  E ) • dx 


=  lA     (J  -  qv).dS 
2(t) 
The  moving  surface  21  (t)  is  to  be  thought  of  purely  mathematically 
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as  a  one-parameter  family  of  surfaces  given  by  equations 
X  =  f(a,[3,t)  with  V  defined  as  the  derivative  dx/dt  at  a  fixed 
point  (a, 3)  of  2_(  t ) .   No  relatlvlstic  or  galllean  conceptions 
are  involved.   Since  v  enters  the  line  integrals  only  in  the 
combinations  v\B«dx  and  vx  E*dx,  we  see  that  any  component  of  v 
tangential  to  P  is  irrelevant.   Similarly,  only  the  component 
normal  to  the  surface  affects  qv«dS.   In  other  words,  the  para- 
raetrization  of  ^  ( or  of  p)  which  defines  a  "fixed  point"  moving 
with  the  surface  is  immaterial. 

The  integral  forms  of  the  energy  and  momentum  equations  (1.5), 
(1,6)  and  (1,8)  taken  on  a  stationary  domain  D  are,  respectively. 


(1,1I|)     ^ 


dt 


(■ 


D 


J=-  B^   +   iXE^)dV  +  ()   i(E  /  B)-dS  =  -   E-J  dV, 


D 


(1.15) 


dt 


^  B^dV  +  &)     i(E/  B)-dS  =  -  J   E-J  dV, 
D   ^        ^S  ^  ^D 


(1,16)     -X^ 


(E/  B)^dV  +  (b  M^^dS ,  = 


(qE  +  JX  B)^dV. 


D 


Jump  conditions  at  a  discontinuity  surface 

Jump  conditions  across  a  discontinuity  surface  can  be 
derived  from  the  various  Integral  forms.   Although  the  same 
results  are  suggested  by  the  differential  equations  they  are 
not  derivable  from  them.   The  integral  formulation  is  basic 
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here,  leading  to  the  differential  eouations  under  assumptions 
of  smoothness  and  to  jiomp  conditions  under  other  assumptions. 
The  precise  natiire  of  the  jump  condition  depends  upon  the  choice 
of  the  variables  assumed  to  be  bOTxnded. 

Let  the  subscripts  n  and  t  designate  respectively  components 
of  a  vector  normal  and  tangential  to  the  discontinuity  surface. 
Prom  (1,11b)  we  obtain,  under  the  assumption  that  B  is  bounded, 

(1.17)  [B^]      =     0, 

i.e.,  that  the  normal  component  of  B  is  continuous  across  a  dis- 
continuity sxxrfsLce.      Prom  (1.11a),  under  the  assumption  that  B 
and  dB/dt  are  bounded,  we  find 

(1.18)  [E^]   =  0, 

namely,  that  the  tangential  component  of  E  is  continuous  across 
a  stationary  discontinuity  sxirface.   Discontinuities  of  E.  may  occur 
at  isolated  Instants,  however.   For  example,  in  a  conducting  fluid 
if  E.  is  assumed  to  be  discontinuous  at  the  Initial  instant,  a 
shock  wave  results  which  immediately  resolves  the  discontinuity 
(of  course,  dB/3t  is  initially  Infinite).   Practically,  an 
Instantaneous  discontinuity  in  E.  may  be  produced  by  the  head-on 
collision  of  two  shocks.   For  a  moving  discontinuity  sxirface  we 
get  a  more  general  condition  from  (1.12), 

(1.19)  [(E  +  v;<B)^]   =  0. 

The  Inhomogeneous  equations  (1.11c)  and  (l.lld)  require  a 
somewhat  different  treatment  for  q  and  J.   If  discontinuities  in 
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E  and  B  are  allowed  then  we  must  allow  infinities  in  q  and  J. 
Specifically,  we  obtain 

(1.20a)     n/  [B^]   =   n  X[B]  =  tij 

(1.20b)       [E^]   =  Q/X, 

where  j  denotes  a  surface  current  and  Q  a  surface  charge  density. 
The  positive  normal  has  been  taken  in  the  direction  from  side 
(0)  to  side  (l)  where  the  jump  in  a  quantity  Z  is  defined  by 
[Z]  =  Z-,_  -  Z^. 

The  energy  and  momentum  equations  do  not  yield  independent 
jump  conditions,  but  these  conditions  are  also  of  interest.  In 
(1.20a)  we  take  the  scalar  product  with  E.  (which  is  continuous 
across  a  stationary  interface)  and  find 

(1.21)     E. j   =   [-  B]  X  E-n. 

The  left  side  in  (1.21)  has  a  plain  interpretation  as  an  energy 
rate  while  the  right  side  is  associated  with  the  Poynting  flux. 

Example  of  the  c\irrent-carrying  loop 

In  illustration,  we  apply  these  formulas  to  the  simple 
case  of  a  single  loop  of  current-carrying  wire  of  finite  thick- 
ness.  For  simplicity,  we  drop  the  displacement  current.   Let  D 
denote  the  entire  domain  consisting  of  the  wire  D,  ,  and  the 
vacuum  domain  Dp,  and  let  S  denote  the  common  boundary  surface 
of  D,  and  Dp.   Using  the  regularity  condition  that  the  total 
electromagnetic  flux  must  be  negligible  in  the  neighborhood  of 
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Infinity,  we  have  the  rate  of  change  of  the  total  magnetic  energy 
from 

and,  for  the  magnetic  energy  in  the  wire  and  vacuum,  respectively, 
<l-23)    M    {     ^B2dV=  (t  lExB-dS 


°2 


(1-2^)    ^/  ^ 


-  B^dV  =-     E«JdV-(piE^B«dS. 
P-  J  J  t^ 

^1  °1  ^ 

Here,  the  positive  normal  to  S  points  into  the  vacuum  domain. 
Insofar  as  E  and  J  can  be  replaced  by  their  average  values  on  a 
cross-section  of  the  wire,  we  may  approximate  the  right  side  in 
(1.22)  by  setting 

(1.25)     j   E»J  dV   =   I  (T)  E'dx 

r 

where  I  denotes  the  total  current  through  a  cross-section  of  the 
wire  (and  is  independent  of  the  choice  of  cross-section)  and 


§ 


E'djc  is  the  e.m.f,  around  the  wire.   This  is  an  electrical  circuit 

C 
approximation.   If  the  wire  is  thin  and  the  current  is  not  too 

strongly  localized  within  a  cross-section,  then  we  may  neglect 

the  magnetic  energy  within  the  wire  compared  to  that  in  the  vacuum. 

We  shall  then  have  in  approximation. 
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I  E.J  dV  =  -  a)  Ex  B-dS   =   10)  E.dx 

D,  '  s  r 


Alternatively,  if  we  assume  that  curl   E  =  0  in  the  wire,  that 
is,  that  B  is  not  changing  within  the  wire,  then  using  div  J  =  0 
we  have  from  MH-V(6.5)'"" 

E«J  dV   =   I  (p  E«dx 

no  matter  what  the  shape  or  size  of  the  wire,   V/e  may  also 
utilize  the  slightly  weaker  assumption  that  E,  is  a  local  surface 
gradient  on  the  surface  S  of  the  wire,  and  obtain  once  more  from 
MH-V(6.19)  that 

(j)-E>B-dS  =  -(j)E«dx   f   ^B-dx  =-I(bE'dx 

s  r        r  r 

where  P  is  a  closed  curve  bounding  a  cross-section  of  the  wire. 

In  the  absence  of  displacement  current,  we  see  that  the 
field  concept  of  Poynting  energy  flux  across  the  surface  of  the 
wire  is  closely  related  to  the  familiar  elementary  circuit  quantity, 
voltage  times  current. 


"Refers  to  Notes  on  Magneto-Hydrodynamics,  N-uraber  V,  equation 
number  (6,5). 
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2,   The  Fluid  Equations 

From  MH-I  (7)  we  take  the  fluid  conservation  equations  in 
Lagranglan  form,  integrating  over  a  domain  D(t)  which  is  carried 
along  by  the  fluid  with  velocity  u. 

(2.1) 


d 

dt 

pdV  =   0 
D(t) 

d 

dt 

pu^dV   +  I 

(2.2)  -^  pu^dV  +  d)  p-^^dS   =     X^dV 

D(t)        S{t)        D(t) 

(2.3)  -^    i      p(S  +  I  u2)dV  +(^  (u^P^^  +  Q^)dS^  =  r  ^dV. 

D(t)  S(t)  D(t) 

The  quantity  X  denotes  the  force  density,  O   the  energy  source,  and 
Q''  heat  flow.   If  we  utilize  Maxwell's  equations  (1.1)  for  the 
electromagnetic  field,  we  shall  specify 

(2.l|.a)    X   =   qE+JxB  +  pg 

(2.1|b)    ^    =     E.J  +  pu-g 

where  g  can  be  interpreted  as  a  gravitational  field  or  any  other 
conservative  force  field.  The  stress  tensor  P"""^  will  usually  be 
assijmed  to  be  scalar, 

(2.5)       P    =  pS'^''^      ^'^^  ^^^   heat  flow  term  is  dropped. 

The  quantity  e  denotes  the  specific  internal  energy  or  internal 

energy  per  unit  mass.   The  equations  (2,1)  -  (2,3)  are  ass-umed  to 
__  .  —  — — 

For  a  relativistic  treatment,  see  Paul  Reichel,  Basic  Notions  of 
Relativistlc  Hydromagnetics,  NYO-7697,  New  York  University, 
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be  supplemented  by  an  eauation  of  state  connecting  the  density  p, 
pressiore  p,  and  Internal  energy  e.   In  these  equations  we  have 
omitted  dissipative  terms  such  as  viscosity  and  heat  conduction. 

The  conservation  equations  (2.1)-(2.3)  may  be  given  in 
differential  form  as  follows: 

(2.6a)     ft  "^  P  ^^^  u  =  0 

(2.6b)    p  ll  +   grad  p   =   qE+JxB  +  pg 

(2.6c)    P  ^  "^   P  div  u   =   (E  +  u/B)«(J  -  qu), 
vjhere  the  Lagranglan  time  derivative  d/dt  is  defined  by 

(2.7)  ^   =  ^  *   (u-V) 

and  all  the  simplifying  assumptions  mentioned  above  have  been 
employed.   The  derivation  of  these  eouations  from  the  integral 
forms  involves  an  amount  of  elementary  manipulation.   An  alternative 
form  of  (2.6c)  in  which  entropy  replaces  internal  energy  as  the 
variable  is  often  useful.   We  use  the  specific  entropy  n     defined 

by 

(2.8)  Td^   =  de  +  pd(l/p) 
to  obtain 

^^'"^  P  If   "   ^(E  +  u/  B).(  J  -  qu). 

In  most  of  the  later  applications  we  shall  not  only  omit  the 
dissipative  terms  of  fluid  dynamics  but  it  will  be  assumed  that 


-  11^  - 


there  is  no  electromagnetic  dissipation  either;  we  only  consider 
it  here  for  purposes  of  illustration.   In  other  words,  v;e  shall 
generally  assiome  the  flow  to  be  adiabatic  and  replace  the  energy 
equation  by  the  condition  that  entropy  is  constant  follo^^fing  a 
particle 

(2.10)     |1  =   0. 

The  total  conservation  equations 

For  a  fixed  domain  and  no  external  force  field  (g  =  0)  the 
equations  (2.2)  and  (2.3)  take  the  form  (cf.  MH-I(6)) 

(2.11a)    ^  I   pu^dV  +  (t)  (pu^u^'  +  p6^^')dS   =  I  (qE  +  J/  B)^dV 
D  S  D 

(2.11b)    ^  f  P(e  +  I  u^)dV  +  0  |p(e  +  I  u^)  +  pj  u-dS 

•J  \^ 


D  S 


J 


r 

E.J  dV 

D 


In  contrast  to  the  mass  conservation  equation,  these  equations 
should  more  properly  be  called  "lack  of  conservation"  equations. 
For  example,  the  rate  of  change  of  momentum  is  given  by  a  surface 
flow  plus  a  distributed  volume  contribution.   However,  on  combin- 
ing the  corresponding  fluid  and  electromagnetic  eauations  we 
should  obtain  a  true  conservation  equation  without  the  volume 
terms.   In  the  case  of  energy,  this  condition  is  satisfied  whether 
we  combine  (2.11b)  with  the  Maxwell  formula  (l.lL).)  or  the  pre-Maxwell 
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formula  (l,l5)  which  omits  displacement  ?iarrent.   For  momentum  we 
obtain  conservation  only  if  we  drop  or  retain  the  force  term 
qE  in  (1.16)  together  with  displacement  current.   Independently 
of  this  criterion,  it  is  easy  to  see  that  in  anything  but  a  very 
poor  conductor,  qE  is  always  small  compared  to  J x  B,   Nevertheless, 
if  displacement  current  is  kept,  then,  whatever  the  physical 
importance  of  the  force  qE,  it  should  be  retained  in  the  equations 
for  mathematical  soundness.   For  the  same  reason  when  displacement 
current  is  dropped  then  the  term  qE  must  also  be  dropped. 

In  summary,  together  with  the  Galilean  Invariant  pre-Maxwell 
equations,  (1,2),  one  should  take  the  force  equation 

(2.12)  P  It  ^   grad  p   =   JXB 

and  either  the  energy  enuation 

(2.13)  p  ^  +   p  div  u  =   (E  +  uyB).J 
or  the  entropy  equation 

A 

(2.11+)  p  |1     =  ^(E  +   ux  B).J. 

The    total    conservation   eouations    are 

(2.15)        ^    I    pu^dV  +  (D   jpu^uJ   +   p5^J'   -  Mg'^ 
D  Q 


dSj     =      0 


I 

,2.16)        ^||p,g.|u2).AB2| 
D 


dV 


+  (^  \p{e   +  ^u^)u^   +   pu^   +  -   (E/B)^/    'dS.   =   0 


S 
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V;ith  the  full  Maxwell  equations  (1.1),  relativistic  fluid  equa- 
tions would  be  more  suitable,  but  as  a  compromise,  one  might  take 
the  force  equation  (2.6b)  and  either  the  energy  equation  (2.6c) 
or  the  entropy  equation  (2,9).   The  total  conservation  equations 
are 


(2.17)    ^ 


and 
(2.18) 


d_ 
dt 


jpu^  +  /(ExB)^r  dV  +  '  Ipu^uJ  +  p6^«5  -  M^^j  dS  .  =  0, 
D  '^  -^  S  ~ 


_d_ 
dt 


J 


pe  +  iu^)  +  ^^   +  ^  E^  i  dV 
2       2ti      2     ; 

+  (p  (p(e  +  |u^)u^  +  pu-^'  +  -(EXB)^  I  dS 


=  0. 


It  is  worth  noting  that  in  the  case  with  no  displacement 
current,  q  occurs  only  in  the  equation  q  =  Xdiv  E;  this  relation 
serves  to  determine  q  after  E  (and  all  other  variables)  have  been 
determined  by  other  conditions.   Conservation  of  charge  will  in 
general  be  violated. 

It  is  customary  to  think  of  the  Poynting  flux  as  the  energy 
passing  through  the  surface  by  radiation  alone,  but  it  actually 
represents  the  entire  flux  of  field  energy.   We  observe,  for 
example,  that  the  appearance  of  the  Poynting  flux  in  the  pre- 
Maxwell  equation  (2,16)  is  not  associated  with  a  transfer  of 
radiated  energy. 

Jump  conditions  at  discontinuity  surfaces 

The  Jump  condition  derived  from  the  mass  conservation  equa- 
tion (2.1)  for  a  stationary  discontinuity  surface  is 
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(2.19)  [pu^]   =   0. 

If  there  is  flow  through  the  interface,  u  7^  0,  we  have  a  shock 

wave.   Shock  discontinuities  are  considered  in  detail  in  MH-VIII. 

For  present  purposes,  we  restrict  ourselves  to  discontinuities 

with  u  =  0  (Helmholtz  discontinuity), 
n  "^ 

From  the  momentum  conservation  eauation  (2.15)  (without 
displacement  current)  and  from  (I.I7)  we  obtain  the  jiomp  condition 

If  B  7^  0  it  follows  that  there  is  no  jump  in  these  variables. 

'  If  B   =  0,  we  have 
n     ' 

(2.20)  tP  +  -^  B^]   =   0»        (B^  =  0)- 

Exactly  the  same  result  holds  for  a  moving  discontinuity  surface. 
This  result  depends  upon  the  assumption  of  scalar  pressures.   In 
this  case,  the  magnetic  force  across  a  discontinuity  surface  is 
balanced  by  a  pressure  force  which  is  normal  to  the  surface.   The 
Maxwell  stress  is  normal  only  when  B  is  either  tangential  or  normal 
to  the  surface,  and  in  the  latter  case,  B  is  continuous.   There 
is  one  important  situation  in  which  the  assumption  of  scalar 
pressure  can  be  dispensed  with.   This  is  the  case  of  a  current- 
carrying  discontinuity  surface  which  separates  a  fluid  region 
where  B  vanishes,  from  a  vacuum  magnetic  field  (cf.  MH-X),   In 
that  case  the  pressure  must  be  exerted  normally  and  satisfies  the 
press-ure  balance  (2,20), 

Prom  the  momentum  equation  (2.17)  with  displacement  current 
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vje  have,    using    (1.18) 

[p]n^   =  i[B^]3^   +  ^[B^ln^  -    X.  eHe]    -   ^[E^Jn^   =   0. 
jj,            n        d\i  t     n  d 

It  follows  on  taking  the  scalar  product  with  n,  that 


(2.21a)    [p  +  -i;^  B^  -  J  E^]  =  0; 


and  taking  the  cross-product  with  n  we  have  from  (1,20) 
(2.21b)    B^( j  /  n)  -  QE^   =   0. 

At  a  perfect  conductor  (see  the  following  section),  we  would 

have  E,  =  0  and  the  relations  (2,21)  are  identical  with  those 

1   2 
without  displacement  current;  namely,  [p  +  p—  B  ]  =  0  and  B  =0, 

The  existence  of  a  surface  current  J  is  incompatible  with  the 

concept  of  finite  conductivity.   Consequently,  for  a  finitely 

conductive  fluid,  B  must  be  continuous  and  we  must  have  either  E 

continuous  or  E,  =  0,   In  the  event  E.  =  0,  then  for  the  problem 

of  determining  the  fields  on  one  side  of  the  discontinuity  surface, 

we  may  as  well  assume  the  conductor  on  the  other  side  to  be 

"perfect." 
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3 

3.   Ohm's  Law; 

If  J  is  an  \inknown  of  a  problem  it  is  necessary  to  supplement 
the  electromagnetic  equations  (1.1)  or  (1.2)  with  an  additional 
vector  equation  which  is  generally  a  counterpart  to  Ohm' s  Law, 
For  a  perfect  conductor,  we  take 

(3.1)  E  +  uy  B  =  0; 

in  this  case  the  fluid  is  non-dissipative  (cf,  (2.9)).   More 
generally,  we  have  (see  MH-IV(12)y, 

(3.2)  E   =  R||  j"   +  CJ  ^  B   +  R,  J-L 

u    II   u    ^  u        1  u 

where  R||  (parallel  resistivity),  R,  (perpendicular  resistivity), 
and  C  (Hall  coefficient)  are  three  phenomenological  "constants'" 
of  the  medium.   The  symbols  E   and  J   denote  values  measured  in 
a  coordinate  system  moving  with  the  fluid  velocity  u,  namely. 


(3.3) 


E   =   E  +  ux  B 
u 


J^   =   J  -  qu  ; 


^^Jhen  displacement  current  is  omitted  we  take  J   =  J.   The  symbols 


u 


II  and  J_  denote  components  perpendicular  and  parallel  to  the  B-field; 

that  is,  J  =  j"  +  J-L  where 
'   u    u    u' 


(3  k)     jll  -  ^^^'^u^  jl  _   3x  (J.B) 

^^'^'  '^n  -     -^2     '       "^u  "     g2 

The  dissipation  (rate  of  entropy  production)  is  defined  by 



See  MH-IV  for  a  general  discussion. 
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(3.^)    E^-J,  =  Rn  (j;;)^  ^    Rx  (Ju)' 


The  complete  set  of  fluid  magnetic  equations  for  a  perfectly- 
conducting  medlTom  (3.1)  form  a  mathematically  tractable  system, 
the  Hall  effect  in  (3.2)  usually  makes  for  complication  even 
though  it  is  non-dissipative . 

Ohm' s  law  is  usually  used  to  eliminate  E  rather  than  J,   In 
particular,  we  have  for  a  perfect  conductor 

(3.6)     ^  -  curKuX  B)  =  0 
dt 

or,  in  integral  form,  (cf.  (1,12)), 


(3-7)     ^ 


\J 


B.dS   =   0. 
(t) 


This  result  means  that  magnetic  flux  is  "frozen"  to  an  element 
of  surface  which  moves  with  the  medium.   Within  a  perfect  conductor 
magnetic  lines  are  identifiable  in  the  sense  that  the  fluid  motion 
carries  magnetic  lines  into  magnetic  lines  (sometimes  it  is  neces- 
sary to  be  careful  in  fixing  magnetic  lines  to  a  conducting  medium; 
compare  the  example  of  MH-V,  p,  129). 

Prom  the  simple  form  of  Ohm's  law,  E  =  RJ  ,  we  obtain,  on 
neglecting  displacement  current, 

(3.8)     —  -  curl(uXB)   =   -  i  curKR  curl  B). 
dt  M- 

In  addition  to  being  carried  with  a  surface  element  the  magnetic 
flux  diffuses  with  a  diffusivity  R/tJ.. 

cf .  K,  0,  Frledrichs  and  H.  Kranzer,  Notes  on  Magneto-Hydrodynamics 
No.  VIII,  NY0-6l|86,  Nonlinear  Wave  Motion. 
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A  simpler  form  of  this  diffusion  equation  is  obtained  in  a 

two-dimensional  geometry  for  which  B  has  a  z-component  only,  u  has 

7 

only  X  and  y  components,  and  no  variables  depend  on  z.    Inter- 
preting B  as  a  scalar,  and  using  (3.2),  we  find 

(3.9)    II  +   div(uB)  =  ^  dlv  (Rj^  grad  B). 

The  parallel  resistivity,  R||  does  not  enter,  nor  does  4  if  it  is 
a  constant. 

When  displacement  current  is  dropped.  Ohm's  law  together  with 
curl  B  =  (iJ  allows  one  to  eliminate  all  electromagnetic  variables 
in  terms  of  B.   In  addition  to  the  fluid  equations  (which  contain 

a  term  curl  B^  B  in  the  momentum  eouation  and  something  like 

2 
(curl  B)   in  the  energy  equation),  there  is  an  eouation  such  as 

(3.6)  or  (3.8),  or  a  slightly  more  complicated  equation  derived 

from  the  general  form  of  Ohm's  law,  (3.2), 


7 

'The  one-dimensional  case  is  analyzed  in  NYO-769I,  Marian  H,  Rose, 

On  the  Diffusion  of  a  Conducting  Fluid  Across  a  Magnetic  Field. 
The  non-dlssipative  problem  is  treated  in  Notes  on  Magneto-Hydro- 
dynamics VII,  Fluid  Dynamical  Analogies  by  A.  A.  Blank  and  Harold 
Grad. 
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h,.      Characteristics 

We  have  considered  the  two  systems  of  non-dissipative  fluid 
magnetic  equations:   a  galilean  system  which  combines  classical 
equations  of  the  fluid  with  the  non-relatlvistic  pre-Maxwell  equa- 
tions, and  a  hybrid  system  which  combines  the  non-relativistic 
fluid  equations  with  the  Lorentz-invariant  field  equations  of 
Maxwell.   In  order  to  evaluate  the  two  systems  we  shall  compare 
them  to  the  complete  relativistic  system  as  a  standard. 

We  have  already  seen  that  for  the  two  systems  under  study  the 
conservation  principles  take  somewhat  different  forms.   More  refined 
distinctions  can  be  made  by  studying  the  characteristics  of  the 
systems,  that  is,  by  examination  of  the  propagation  of  weak  dis- 
tTorbances,   These  systems  and  the  relativistic  system  are  all 
symmetric  hyperbolic.^   The  existence  theory  and  general  properties 
of  such  systems  have  been  extensively  developed. 

In  each  of  these  systems  the 


k 
For  the  galilean  system,  see  K.  0,  Priedrichs,  Nonlinear  Wave 

Motion  in  Magneto-hydrodynamics,  LAMS-210^.   K.  0,  Priedrichs  and 

H,  Kranzer,  Nonlinear  Wave  Motion,  Notes  on  Magnetohydro dynamics  VIII, 

NYO-61^.86,  New  York  University,  Institute  of  Mathematical  Sciences, 

(1958). 

For  the  relativistic  system,  see  P,  Reichel,  Basic  Notions  of 
Relativistic  Hydromagnetics,  NYO-7697,  New  York  University,  Institute 
of  Mathematical  Sciences,  (1958). 

Also,  B.  Zumino,  Some  questions  In  relativistic  hydromagnetics, 
Phys.  Rev,  108,  Noo  5,  1116-1121  (1957). 
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characteristic  equation,  which  has  as  Its  roots  the  characteristic 
velocities  at  which  weak  signals  propagate,  is  a  polynomial  of 
sixth  degree.   Letting  V  denote  the  characteristic  velocity  as 
measured  in  a  frame  at  rest  with  respect  to  the  fluid  we  obtain 
for  the  characteristic  equation  in  the  non-relativistic  case, 

(l+.l)    [^pV^  -  B^lLix^p^V^  -  iapV^(tipa^  +  B^)  +  tipa^B^]  =  0 

and  in  the  relativlstic  case, 

^)  -  B2][vVp-^p-==-  -.^^ 
0  c 


(1^.2)    [V^ltip'"'"  +  %)  -  b2][vVp'''(p'''  +  %)  -  V^  iy.p''{^pa.^   +  B^ ) 

2  2 

a'^B  1      p  p 

+  tip  —^     +  [ipaB^    ]  =  0, 
c   J 


2        -1 

In  these  equations  n  denotes  the  wave  normal,  c   =  ( M-  k)         is  the 

squared  light  speed  and   p"  =  p  +  —^    (p  +  pe)  where  p  is  the  rest- 

0 

mass  density  (rest-mass  is  conserved),  and  e  represents  the  internal 
energy  per  \anlt  rest  mass  measured  in  the  frame  at  rest  with  respect 
to  the  fluid. 

The  principal  difficulty  with  the  non-relativistic  system  is 
that  in  the  low  density  range,  some  of  the  characteristic  speeds 
become  unbounded.   In  fact,  for  the  roots  of  (L|-,l)  we  have 
asymptotically  in  p: 

P   a^B^    B^   „2 

(l|.3)      V2 ^    n   B   . 

g2   '  tip  *  tip 

In  other  words,  for  sufficiently  low  densitl-es  the  non-relativistic 
equations  permit  the  propagation  of  a  signal  in  excess  of  light 
speed.   For  comparison  we  observe  that  the  corresponding  roots  of 
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(l4-,2)  are  given  by 

, ,  ,  s    ,,2      n       n p 

(1^.1^)    V  ~  — 2~»  2TT  '   27~2  • 

B'^  ixp  +  BVc     tip  +  B  /c^ 


Prom  these  observations  it  is  clear  in  what  realm  the  non-relativistic 

system  is  an  adequate  approximation  to  the  relativistic  system. 

At  this  point,  it  is  reasonable  to  ask  whether  the  hybrid 

system  has  any  advantages  to  compensate  for  the  lack  of  invariance. 

In  particular,  to  be  useful,  the  hybrid  system  ought  to  be  simpler 

than  the  relativistic  system  and  at  the  same  time  provide  a  better 

approximation  than  the  non-relativistic  system  in  the  low  density 

limit.   For  the  hybrid  system,  the  propagation  speeds  are  taken 

with  respect  to  a  stationary  frame  since  the  concept  of  velocity 

relative  to  the  fluid  has  no  invariant  meaning.   We  use  V   to 

designate  velocity  relative  to  the  given  frame  and,  for  comparison 

with  the  other  systems,  set  V  =  V   -  u  .   The  characteristic  equation 

•^  '  on 

for   the  hybrid   system  takes   the   rather   complicated   form 

where 

f  =   txpV^    -    B^,      ^=   t,2pV(v2-a2)    -    /I,    n  =   ixpCV^B^-aV) 

B^  =  B^  -  B^  ,   U2  =  U'(B-nB^)/B^,   u^  =  u-(nx  B)/b^. 
The  product  J  '^  is  the  left  side  of  equation  (l^..l).   It  is  clear 
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that  this  is  a  far  more  complicated  equation  than  either  of  (I4..I) 
and  (11,2).   In  the  zero  density  limit  the  characteristic  speeds  are 
given  by 


(1+.6) 


r  (V  -  u  )^  =  V^  =  a^  ^  , 
on  g2  ' 


< 


B 


V^  -  u   =  V 
o    n 


=  -w  -1.  ^/ 


c  -  u 


^1 


V^  =  (V  +  u  )2  =  c^ 
o         n 


where  u 


II 


U'B 
iBl 


and  u' 


■1 


2    2 
u  -  u,, 


It  is  clear  that  the  propagation  speeds  (in  either  the  rest  frame 
or  the  fluid  frame)  are  not  bounded  by  light  speed  c  unless  the  fluid 
velocity  is  dropped  as  being  small  compared  to  c.   The  hybrid  system 
is  useful  only  when  l/c   is  considered  small  as  well  as  p  and,  with 
that  understanding,  only  the  lowest  order  terms  are  pertinent.   The 
relativistic  correction  to  the  galllean  system  provided  by  the 
hybrid  system  is  applicable  only  in  a  limited  range  and  at  the  cost 
of  greater  complexity  than  even  the  fully  relativistic  system. 
There  would  then  seem  to  be  no  practical  value  in  the  use  of  the 
hybrid  system. 
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